
Fig. 1 : Schematic picture of the surface layer with length scale a 
and Dirac cone of the surface layer 




Fig. 2: Sonderheimer oscillation (SO) in magnetoresistance 




H(T) 

Fig. 3: Dependence of SO on the Fermi energy and tlie disorder strength 




Fig. 4: MR and Hall resistance of sample #3 and #4 
with theoretical fitting curves 
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Abstract 

Topological states of matter challenge the paradigm of symmetry breaking, characterized by 
gapless boundary modes and protected by the topological property of the ground state. Recently, 
angle-resolved photoemission spectroscopy (ARPES) has revealed that semiconductors of Bi2Se3 
and Bi2Te3 belong to such a class of materials. Here, we present undisputable evidence for the 
existence of gapless surface Dirac fermions from transport in Bi2Te3. We observe Sondheimer 
oscillation in magnetoresistance (MR). This oscillation originates from the quantization of motion 
due to the confinement of electrons within the surface layer. Based on Sondheimer's transport 
theory, we determine the thickness of the surface state from the oscillation data. In addition, we 
uncover the topological nature of the surface state, fitting consistently both the non-oscillatory 
part of MR and the Hall resistance. The side-jump contribution turns out to dominate around 1 
T in Hall resistance while the Berry-curvature effect dominates in 3 T ~ 4 T. 

PACS numbers: 
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Symmetry breaking is the paradigm in not only classifying quantum states of matter but 
also describing phase transitions between them, where the correlation length of fluctuations 
between local order parameters diverges at the critical point of a continuous transition 
l|, |2|. On the other hand, topological states of matter are classified by topological quantum 
numbers js, 4|, associated with gapless boundary electronic states and are protected from 
the topological properties of the ground state |5|, l6| . Instead of a divergence in correlation 
length, topological phase transitions are accompanied by changes of the gapless boundary 
modes. One possible mechanism for this phenomenon is the length scale of the boundary 
state becomes of the same order as the bulk size, causing the gapless modes in opposite 
boundaries to be mixed and making such boundary modes gapped In this case, the 
length scale for the boundary mode plays basically the same role as the correlation length 
for the phase transition, which is the fundamental length scale for a topological phase (Fig. 
la). 

Recently, the semiconductor material Bi2Se3 and Bi2Te3 have been verified to be three 
dimensional topological insulators. These insulators are regarded as a novel quantum state 



of matter 



where gapless surface electrons are uncovered in ARPES 12l-ll4|. Although 



the surface states in these materials and the electronic structure in graphene are described by 
Dirac theory (Fig. lb), the surface state of topological insulators is profoundly different from 



the electronic structure of graphene 



15| . This difference originates from the absence of both 



the =ublattice =yn„net.y and valley degenera^. The dhect.on of =pm i= locked with that ot 

the momentum in surface Dirac electrons This completely suppresses backscattering 
due to time- reversal invariant impurities, allowing a super-metallic state |l^. 

In this letter we focus on MR and Hall measurements, both of which are of high impor- 
tance for the fundamental understanding and practical applications of topological insulators. 
We observe an oscillatory behavior in MR at low magnetic fields of up to 4 T. This behaviour 
can be identified as Sondheimer oscillation 18(|, where the oscillation period is linearly pro- 
portional to the magnetic field (Fig. 2). Sondheimer's transport theory 18|, combined with 
Dirac dispersion, enables us to determine the fundamental length scale from our experi- 
mental data, which turns out to be about 5 atomic layers. In this respect the Sondheimer 
oscillation can be regarded as an inevitable result of the surface state. The nature of the 
surface state can be explained by the single Dirac-fermion theory. This explains not only 
the non-oscillatory part of MR but also the topological properties implicit in Hall resistance 
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in a quantitative and consistent way (Fig. 4). In particular, we show that the Hall resis- 
tance of the surface state is dominated by the side-jump around 1 T and below while the 
Berry- curvature contribution is dominant at higher fields (Fig. 5 in SI III). 

In our measurements we used defect-controlled Bi2Te3 single crystals. Usually, as-grown 
Bi2Te3 single crystals are p-doped because of the anti-site defects in Bi sites |27|. In order to 
tune the Fermi level, we have controlled the amount of defects by adding extra Bi or Te; the 
doped Bi tends to increase the anti-site defects, while the doped Te tends to decrease them. 
Based on this strategy, we have succeeded in growing a range of Bi2Te3 single crystals, from 
fully p-doped to fully n-doped regions. Though rare, we have obtained almost insulating 
Bi2Te3 single crystals. The chance for obtaining these samples is around 3 %. The carrier 
type is determined by thermoelectric power at room temperature and also by Hall sign 
measured at 4.2 K. Samples jjl and {12 are hole-doped while sample jjS is electron-doped. 
Samples jjS and tj4 are close to insulators, as such, they are expected to show the topological 
properties of the surface state well. See the supplementary material (SI I) for our sample 
preparation. 

Magnetoresistance (MR) and Hall effect measurements have been carried out by a six- 
probe method at 4.2 K using a superconducting magnet up to 4 T and a 60 T pulse magnet 
at ISSP in Tokyo university up to 55 T. Here, the direction of the magnetic fields is applied 
perpendicular to the naturally cleaved plane, on which the current is applied. For the 
MR and Hall measurements, we carefully contacted the lead wires reduce the induction 
noise for high-field pulse-magnet experiments. We have taken the ant i- symmetrized and the 
symmetrized parts as Hall and longitudinal resistances, respectively. 

ARPES has unveiled only a sing 
Therefore, in order to analyze our experimental data, we introduce an electromagnetic vec- 
tor potential A and a Zeeman term into the single Dirac-fermion theory. In addition, we 
take into account impurity scattering at the level of a Born approximation. Our theoretical 
analysis reveals that the orbital contribution or the effect of the vector potential on exper- 
imental data is irrelevant in the region of magnetic fields below 4 T (SI IV). However, its 
influence on MR and Hall resistance can appear at higher magnetic fields, where Landau 
levels are fully developed. 

We model the surface state as a thin layer with thickness a (Fig. la), which is used in the 
Boltzmann equation (SI II). This approach is essentially the same as what Sondheimer used 
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with metallic thin films except for the band structure, where non-relativistic electrons are 
replaced with Dirac fermions in the presence of the Zeeman term. The main consequence of 
this is the quantization of motion along the direction normal to the surface which produces 
an oscillatory component of MR. In our measurements at magnetic fields below 4 T, the 
periodicity observed is linearly proportional to H [18]. This oscillation is distinguished from 
Shubnikov-de Haas oscillation due to the formation of Landau levels, where the periodicity 
is proportional to 1/H. See Fig. 2. 



Sondheimer's transport theory [18| with the Dirac dispersion gives the following expres- 
sion for the oscillating part of the longitudinal resistivity 

P{H, T) 1 , 



p(0,0) K 



ms). (1) 



where the kernel 



with s = K + iP results from the distribution in the hard-wall boundary condition for the 
2;-direction (SI II). Two parameters appear in this expression, n = a/l and (3 = a/rc, where I 
is the mean free path and is the magnetic length (proportional to 1/H). a is the thickness 
of the surface state. This transport theory produces the if-linear periodicity. It is worth 
noting that the periodicity in MR depends only on f3 while k modifies the amplitude of 
oscillation as shown in Fig. 3. 

Our experimental data shows that the Sondheimer oscillation turns into the Shubnikov-de 
Haas oscillation above 3 T ~ 4 T (Fig. 2). This is consistent with several recent transport 



measurements that show Shubnikov-de Hass oscillation beginning at around 4 T [19|-|22|. 
The appearance of the Shubnikov-de Haas oscillation is the origin for the mismatch of the 
oscillation amplitude. 

We point out that /3 can be written as /3 = ^{kfa)h{kf/k) with the surface thickness a, 
where h = huj^/Ep is the ratio of the Zeeman energy Tilol to the Fermi energy Ep at the 
surface, and kf and k are the Fermi wave- vector and the average momentum, respectively. 
See the supplementary information for details. An important aspect is that h also enters 
both the non-oscillatory part of MR and Hall resistance. Therefore, the actual value of 
h influences not only the periodicity of the Sondheimer oscillation but also both the non- 
oscillatory part of MR and Hall resistance. Precisely speaking, the surface thickness cannot 



be determined from the Sondheimer oscillation alone. Combined with the longitudinal and 
transverse resistances, we can optimize the thickness a and the coefficient 7 simultaneously 
(SI II), where 7 is the ratio TujJl/ Ep at = 1 T. It is interesting to note that the optimized 
7 almost coincides with the bulk value. Furthermore, this 7 value seems to be common to 
both Bi2Te3 and Bi2Se3 at a given field (SI II). Our fitting for the oscillatory part of MR, 
performed consistently for both the non-oscillating part of MR and Hall resistance, gives a 
result for the surface thickness of approximately 5 atomic layers. This is quite remarkable 



in that this value is consistent with that in MBE-grown Bi2Te3 thin films |23 |. 

Next, we focus on the topological nature of the surface state. Strictly, the role of the single 
Dirac-fermion theory is not essential in Sondheimer oscillation although Dirac dispersion is 
utilized. It might be the case that the surface state is realized due to the good surface 
quality of our samples. However, we will show that the Sondheimer oscillation is a signature 
of surface Dirac electrons in Bi2Te3, verifying that the Hall resistance originates from the 
anomalous Hall effect of Dirac theory. In addition, we show that the side-jump contribution 
dominates at low magnetic fields and the Berry- curvature effect dominates at high magnetic 
fields. 

The single Dirac-fermion theory gives an analytic expression for the longitudinal conduc- 
tance 24 1 



,„(^.^^0)^.^^, (3) 

where a = ^"'^ J-^ is a dimensionless parameter which measures the strength of disorder 

'm*nj[Vj ] kf 

with an impurity density nj and an impurity potential Vj^\ while f/ is the Dirac-fermion 
velocity with the Fermi momentum kf. h is the dimensionless magnetic field, introduced in 
Sondheimer oscillation. 



The same Dirac theory results in Hall conductance 24 1 



a,y{H, T^O) = ag'{H, T ^ 0) + a^^{H, T ^ 0), 
at{H, T ^ 0) = a^(H, T ^ 0) + af,^(//, T ^ 0) 



+ aff(iJ,T^O). (4) 

Hall conductance consists of two contributions. The first results from electrons near the 
Fermi surface, referred as normal Hall conductance while the second contribution comes 
from a fully occupied band, called anomalous Hall conductance. We use the normal Hall 



conductance from the Boltzman equation approach for the Sondheimer oscillation, which is 
given by 



The anomalous Hall conductance is also composed of two contributions. The first comes 
purely from the topological character of the band structure, identified with the Berry- 
curvature term a^y{H,T 0) — — ^yj=p, while the second originates from scattering 
with disorder in the presence of the spin-orbit interaction. This disorder contribution is 

separated into the side-jump term a^^{H,T -)■ 0) = -^-^^^j^ + jp^^^'^ and the 
skew scattering term a^^{H,T — )■ 0) = (i+ih^)^ ■ interesting to observe that the 

side-jump term does not depend on the disorder strength. The dimensionless parameter 



= ] y^f in the skew scattering term measures the disorder strength in the third order, 
where Vj^^ is a disorder strength of the third order. 

Based on Eqs. (3) and (4), we obtain the longitudinal and Hall resistances as follows 



[(7,,(i/,T^0)]2 + [a,,(i/,T^0)]2' 



p^y{H, T ^ 0) - ^^^^^^^ ^ ^^J^^ ^ ^ ^^^^ . (5) 

It should be noted that a^y in the denominator cannot be ignored in this case because this 
term is comparable to axx- Here we have two dimensionless parameters, a and rj. However, 
the contribution from the skew scattering turns out to be neghgible (SI HI). Only one fitting 
parameter a remains for both MR and Hall resistance. 

When either holes (sample p and jj2) or electrons (sample jlS) are heavily doped, MR 
curves greatly deviate from the single Dirac-fermion theory (SI III). On the other hand, 
nearly insulating samples (sample i|3 and }|4) display reasonable matches between experiment 
and theory, see Fig 4. The Hall resistance also shows deviation from the single Dirac-fermion 
theory for heavily doped samples (SI III) but not much for nearly insulating samples (sample 
US and ii4), implying that insulating samples are explained by the theory in a quantitative 
and consistent way. These results provide a compelling evidence for surface Dirac electrons. 
In particular, the dominant contribution in the Hall resistance turns out to be the side-jump 
mechanism at fields below 1 T and the Berry-curvature effect at higher fields. See Fig. 5 in 
SI HI. 
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It is also worth noting that the curvature of the Hall resistance in sample {14 is larger 
than that in sample {13. According to Dirac theory with disorder, two parameters affect the 
shape of the Hall resistance: the disorder strength a and the parameter 7 that measures the 
distance from the Dirac point. In Fig. 6 of SI HI, we show how these parameters influence 
the curvature of Hall resistance. By decreasing 7, the Hall resistance becomes straighter 
because the anomalous Hall effect weakens. 

In this letter we have measured the fundamental length scale of the topological insulator, 
the thickness of the surface state, from the Sondheimer oscillation in magnetoresistance. 
This surface state is described by the single Dirac-fermion theory. The topological nature of 
this is verified by the fact that the Hall resistance mainly results from the anomalous Hall 
effect of Dirac theory, which in turn is dominated by both the side-jump mechanism and 
the Berry-curvature effect. 

The surface thickness will diverge at the critical point of a phase transition from a band 
insulator to a topological insulator. Such a phase transition was demonstrated in the HgTe 



quantum well structure when the size of the quantum well was tuned [25|, |26l]. On the 
other hand, the topological phase transition has not yet been achieved in three-dimensional 
topological insulators. Our measurement for the surface thickness can be utilized as an 
important tool, revealing the mechanism of such a topological phase transition. 
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FIG. 1: A schematic picture of the surface layer with a length scale a and the Dirac cone 

at the surface state, a. The schematic diagram for Bi2Te3 shows the surface layer with thickness 
a. The surface thickness a is determined from Sondheimer oscillation in magnetoresistance, orig- 
inating from the quantization of motion within the surface layer, b. The Dirac dispersion of the 
surface state gives rise to topologically nontrivial physical properties. In particular, the anomalous 
Hall effect due to the Dirac cone turns out to dominate in the Hall resistance. See Fig. 5 in the 
supplementary information. 

FIG. 2: Sondheimer oscillation in magnetoresistance. a. The second derivative of the 
resistance with respect to the applied magnetic field shows oscillation with a periodicity in H, 
compared to the theoretical curve (red thick line) based on Sondheimer's transport theory. The 
experimental periodicity matches well with the theoretical curve, but it deviates from theoretical 
values around H = AT. h. Peak and dip number vs. peak and dip position (magnetic fields) in (a). 
Data points arc located on a straight line at low magnetic fields, confirming the H linear periodicity 
instead of the 1/H periodicity. On the other hand, the experimental data will deviate from the 
straight line in high magnetic fields due to the appearance of Shubnikov-de Haas oscillation. The 
second derivative of the magnetoresistance measured up to 55 T is plotted with respect to H (c) 
and 1/H (d). This comparison reveals that the Sondheimer oscillation exists in the region of low 
magnetic fields while the Shubnikov-de Haas oscillation with a periodicity in appears at high 
magnetic fields. 

FIG. 3: Dependence of Sondheimer oscillation on the Fermi energy and the disorder 
strength, a. Dependence of Sondheimer oscillation on the Fermi energy. ^ = h/H measures 
the distance of the Fermi energy from the Dirac point, where h = hu)L/Ep is the dimensionless 
magnetic field given by the ratio of the Zeeman energy Hul to the Fermi energy Ep. Increasing 7, 
i.e., as the Fermi surface becomes close to the Dirac point, the period of the Sondheimer oscillation 
decreases, b. Dependence of Sondheimer oscillation on the disorder strength, k = a/l measures 
the mean free path. It does not affect the periodicity, changing the amplitude of the oscillation 
only. c. Peak and dip number vs. peak and dip position (magnetic fields) as a function of 7 with 
a fixed K. This confirms our conclusion in Fig. 3a. 
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FIG. 4: Magnetoresistance and Hall resistance for sample [j3 and with theoretical 
fitting. The magnetoresistance and Hall resistance of sample (jS are displayed in a and b, respec- 
tively, together with theoretical curves (red thick line) based on the single Dirac-fermion theory. 
The same quantities of sample p are presented in c and d. We emphasize that the theoretical 
curves for Hall data are based on the parameters from our fitting of magnetoresistance data. In 
other words, 7 (the ratio of the Zeeman energy to the Fermi energy) and a (the disorder strength 
in the longitudinal resistance) are determined from both the Sondhcimer oscillation and magne- 
toresistance completely, and there are no free parameters for the Hall resistance. See the text. 
Reasonable matches between the experimental data, particulary for the Hall resistance, and the 
theory reveal that the nature of the surface state is described by the single Dirac-fermion theory. 
Furthermore, we found that the Berry-curvature term dominates the experimental data of Hall 
resistance around 3 T ~ 4 T while the side-jump mechanism works around H = 1 T, confirming 
the topological origin of the transport phenomena in the surface state of topological insulators. 
See Fig. 5 in the supplementary information. 



11 



Supplementary information for Sondheimer oscillation as a 
fingerprint of surface Dirac fermions 

Heon-Jung Kim^'*, Ki-Seok Kim^'^, Mun Dae Kim"^, S.-J. Lee^, J.-W. Han\ 
A. Ohnishi^, M. Kitaura^, and M. Sasaki^'^ A. Kondo*^, and K. Kindo*^ 
^Department of Physics, College of Natural Science, 

o ■ 

CN ■ Daegu University, Gyeongbuk 712-714 Republic of Korea 

, '^Asia Pacific Center for Theoretical Physics, 

POSTECH, Pohang, Gyeongbuk 790-784, Korea 

^Department of Physics, POSTECH, 

D ■ Pohang, Cyeongbuk 790-784, Korea 

! ^Institute of Physics and Applied Physics, 

\ Yonsei University, Seoul 120-749, Korea 

I ^Department of Physics, Faculty of Science, Yamagata University, 

^ ■ Kojirakawa 1-4-12 Yamagata, 990-8560, Japan 

O ^ 

, ^Institute for Solid State Physics, University of Tokyo, 

Kashiwanoha 5-1-5, Kashiwa, Chiba 277-8581 Japan 
(Dated: January 20, 2013) 

Abstract 

In this supplementary information we discuss how to fit our experimental data to theoretical 
predictions of transport on the surface of Bi2Te3. First, we describe the sample synthesis and its 
resulting character. Second, we present the detailed procedure for fitting Sondheimer oscillation 
data to theory, where the fundamental length scale of the surface state is determined. Third, we 
discuss the non-oscillatory part of magnetoresistance and Hall resistance based on the single Dirac- 
fermion theory, where the topological nature of the surface state is verified. Last, we consider the 
quantum Hall effect by examining electromagnetic vector potential. The theoretical result from 
the quantum Hall effect shows that our experimental regime is far from the quantum Hall effect 
and thus, the effect of the electromagnetic vector potential is negligible in our analysis. 

PACS numbers: 
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I. SAMPLE PREPARATION AND CHARACTERIZATION 



Single crystals of Bi2Te3 were grown by a modified bridgeman method, where Bi2Te3 
powder is melted and crystalized in an evacuated quartz ampoule several times by slow 
cooling. The sample was cooled from 850° C to 550° C with a cooling rate of -10 K/h. Since 
as-grown Bi2Te3 is usually p-doped (p-Bi2Te3) due to the anti-site defects at the Bi sites[l|, 
we have controlled the amounts of these defects by adding extra Bi or Te; the doped Bi tends 
to increase such defects, while the doped Te tends to reduce them. To grow topologically 
insulating samples (TI-Bi2Te3), we should add small amounts of Te to the p-Bi2Te3 ingots. 
In our experiments, we found that 7 ~ 8 % extra Te was necessary in order to synthesize 
TI-Bi2Te3. Even if such an ingot were carefully grown, only few single crystals from that 
ingot showed the desired nonmetallic T-dependent resistivity curves, indicating that large 
ingots usually contain weak concentration- or defect-gradient Q]. Therefore, we had to select 
and check every crystal from the large ingots to confirm their non-metallic characteristics. 
In the selection process, we measured the temperature-dependent resistivity down to 77 K 
and thermoelectric power at room temperature of about 200 samples. The n-Bi2Te3 single 
crystals were obtained by the further doping of Te to produce the TI-Bi2Te3 ingots. 

For our experiments, we selected two p-Bi2Te3 samples (fully p-type jjl and lightly p-type 
jj2), two TI-Bi2Te3 samples ({13 and tj4), and one n-Bi2Te3 sample (tl5). The carrier sign was 
determined from the sign of the thermopower at room temperature and was also checked 
using Hall resistivity at 4.2 K, both of which gave the same sign. Magnetoresistance (MR) 
and Hall effect measurements have been carried out by a six-probe method at 4.2 K, using 
a superconducting magnet up to 4 T and a 60 T pulse magnet at ISSP in the university of 
Tokyo up to 55 T. Here, the direction of the magnetic field applied is perpendicular to the 
naturally cleaved plane, on which the current is applied. For MR and Hall measurements, 
we carefully contacted the lead wires to reduce induction noise for high-field pulse-magnet 
experiments. We have taken the anti-symmetrized and symmetrized parts in Hall and lon- 
gitudinal resistances, respectively. 

Figure 1 shows the temperature dependency of the resistivity for the p-, TI-, and n-Bi2Te3 
single crystals. The resistivity for the both p- and n-type samples (jjl, tl2, and tl5) decreases 
monotonically with decreasing temperature, which is a typical metallic characteristic, while 
the resistivity for the TI samples (jjS and tl4) increases below ~180 K and then tends to 



2 



25 



#1 [x5] 




T {K) 

FIG. 1: Temperature dependency of resistivity for samples fjl - jjS. 

saturate below 50 K. The magnitude of the resistivity for the TI-Bi2Te3 single crystals is 
5 ~ 10 times larger than those for the p- and n-type samples. The non-metallic nature 
observed for the TI samples is consistent with those reported by the Princeton group j2|. 
The oscillating part of MR was separated by fitting the background with a polynomial. 

We have analyzed all five samples and discussed them in both the main text and the 
supplementary information. Among them, the TI samples ({jS and tl4) are best fitted to 
the single Dirac-fermion theory while the fitting becomes worse for the p- and n-Bi2Te3 
samples. Deviation between the data of samples jjl, tl2, and jjS and the theoretical curves 
can be seen in Fig. 3 and Fig. 4. These figures clearly show that the linear part in MR of 
the experimental data deviates from the theoretical fitting significantly, implying that bulk 
channels for conduction are important beyond the single Dirac-fermion theory. Judging from 
our analysis, the Fermi level of the surface states for sample tl3 and {14 lies below the Dirac 
point, corresponding to hole doping, while that in the bulk is between the bulk gap. The 
difference of the Fermi level between the surface and bulk originates from a band bending 
as reported before jsl, |^ . On the other hand, the bulk Fermi level is below the valence band 
maximum for sample [jl and {12. This explains the metallic properties and the disagreement 
between theory and experimental data for samples jjl, tl2, and jjS. 
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II. SONDHEIMER OSCILLATION 



Sondheimer oscillation appears in metallic thin films, where a finite thickness of a surface 
state gives rise to a quantization of motion, resulting in an oscillation in MR and Hall 
resistance jsl- While Shubnikov-de Haas oscillation displays periodicity proportional to 
the periodicty of Sondheimer oscillation is proportional to H. This periodicity allows us to 
find the thickness of the surface state. Although Bi2Te3 is a bulk material (semiconductor), 
the topological structure of the ground state gives rise to a surface state with thickness 
a, which is protected from time-reversal invariant perturbations 0-9]. In this respect, the 
Sondheimer oscillation is undisputable evidence for existence of a surface state. 

We start from the Boltzmann equation 

-l(E + -xH)-Vuf + v- Vrf = -l^. (1) 



^ V c 

/ = fo + fi{v, z) is a non-equilibrium distribution function with its equilibrium part /q. The 
non-equilibrium part /i depends on the z coordinate, v = hk/m* is the average velocity, 
where m* is an effective mass of the surface Dirac electrons, and k is the average momentum, 
which is determined later. The dispersion is given by = —fi + \Jv'jk'^ + {g*Hy in the 
presence of a z-directional magnetic field H, where f/ is the Dirac velocity and the Fermi 
momentum is kf = fi is the chemical potential at the surface, and g* is the 

Lande-g factor of the surface electron. H = Hz is an applied magnetic field in the z 
direction and E = E^x + Ey-y is an electric field, where a y-directional electric field is 
induced, r is the mean-free time, which measures the strength of disorder. 



Following the same procedure as in the original paper of Sondheimer 



we obtain 



p{H,T) = P^ms), (2) 

where p{H, T) is the resistivity. 

Po — is the residual resistivity with the density n = 
'^^ surface electrons, and k = j, where I = vt is the mean free path and a is the 
surface thickness, determined from fitting. results from the distribution function in the 
hard-wall boundary condition for the z direction, given by 

1 _ 1 _ 3 ^ r , -su( 1 _ M (o^ 
<P{s) ~ s 8^2 + 2.2 7, W U^J 

after integrating over z. s = n + if3, where [i = f- with the magnetic length = 



For the numerical analysis, it is important to express the ^ variable in terms of dimen- 
sionless parameters because the surface thickness a is determined from an appropriate choice 
of one parameter in ^, referred to as 7, which will be discussed later. It is given by 

1 .hwikf 



where Hujl — ^1^^ is the effective Zeeman energy and Ep — is the Fermi energy with 



1.2 

"2m*c ^" ""^ — — - -"--6J -i- — 2m* 

an effective mass m*. k/kf is determined from 



I ,,2_ / 

k 1 t c I — 

where we introduce several dimcnsionless parameters, scaled by the Fermi energy, such that 
an effective dispersion e = vk/Ep with v — the Zeeman energy h — fuJi/EF, the 

chemical potential n' — ji/Ep — vT+7?, and an effective temperature t — T/Ep. 

Based on this formulation, we fit the oscillation data of MR. See Fig. 2 for sample ji2. 
First, we determine k ^ 0.02, which characterizes the strength of the disorder for the best 
match of the amplitude with the Sondheimer oscillation. It is of note the matching cannot 
be perfect. In fact, the variation of k does not modify the periodicity of the Sondheimer 
oscillation although it changes the oscillation amplitude. 

The surface thickness is determined from the periodicity AH in the experimental data. 
Considering that Eq. (5) fixes the ratio of k/kf, two parameters in (3 remain unknown, 
which correspond to kfa and 7 which is the ratio fkjj/ Ep at = 1 T. Therefore, we cannot 
determine the thickness of the surface state from the oscillation period alone. As discussed 
in the main text, the same value of 7 is utilized for both the non-oscillatory part of MR and 
Hall resistance. If the value of 7 were incorrectly chosen, we could not fit both the non- 
oscillatory part of MR and Hall resistance consistently. See Fig. 6-b. Combined with both 
longitudinal and transverse resistances, we found that 7 ~ 0.44 fits the experimental data 
well. This value is quite interesting because it almost coincides with the 7 value determined 
from the values of the bulk parameters. For example, if we use bulk values of the Lande-^f 
factor and effective mass, which are g* ^ 13.7 and m* ~ O.lmg (mg is the bare mass of an 
electron), we obtain Itlul ~ 100 K at = 1 T, giving rise to /i ~ 0.44 with a typical bulk 
value for the Fermi energy Ep 230 K. Since we cannot determine g*, m*, and Ep at the 
surface from our experiment, each value at the surface may not be the same as the bulk one. 
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FIG. 2: a. Sondheimer oscillation for sample p. The periodicity is well matched, but the oscillation 
amplitude deviates rather a lot due to the appearance of the Shubnikov-de Haas oscillation above 
4 T. Sondheimer oscillations are also observed for the p- and n-doped samples, where the surface 
channel for conduction coexists with the bulk conduction, b. Peak and dip number vs. peak and 
dip position (magnetic fields) in (a). Data points are located on a straight line, confirming the H 
linear periodicity instead of the 1/H periodicity. 

However, it should be noted that our 7 seems to be universal in both Bi2Te3 and BioSes 
although both g* and Ep in Bi2Se3 are almost four times larger than those in Bi2Te3 |3|. 

III. MAGNETORESISTANCE AND HALL EFFECT 

Theoretically, the longitudinal resistance results from the transport of electrons near the 
Fermi surface, and is described quasi-classically or quantum-mechanically. In this case, the 
quasi-classical treatment based on the Boltzmann equation gives the same result as the 
quantum mechanical treatment based on the Kubo-formula. On the other hand, there are 
various contributions with a topological origin in the Hall resistance, which is beyond the 
conventional treatment used in the Boltzmann equation approach. The Boltzmann equation 
needs additional terms in order to mimic the Kubo-formula [121, 13] • 



Since several terms exist in the Hall resistance, which are explained in the main text, 
an important question is which terms give a large contribution to the Hall resistance. We 
explicitly demonstrate that the side-jump mechanism dominates in low magnetic fields while 
the Berry-curvature effect dominates the Hall resistance at high magnetic fields (Fig. 5). In 
contrast, the skew scattering term turns out to be negligible. We also check the dependence 





100 




FIG. 3: Magnetoresistance for a. ftl, b. jJ2, and c. [15, where the disorder strength is utihzed as a 
fitting parameter, given by q ~ 0.4, a ~ 0.5, and a ~ 1.5, respectively. As discussed in section I, 
the presence of bulk conduction channels does not allow us to describe these samples purely within 
the single Dirac-fermion theory. 

of the Hall resistance on either the disorder strength a (Fig. 6-a) or the parameter 7, which 
determines the Fermi level (Fig. 6-b). Roughly speaking, to vary 7 is to rescale the x-axis 
while to change a is to rescale the ?/-axis. Figure 6-a and Fig. 6-b show how these parameters 
affect the shape of the Hall resistance. These results confirm that if we do not choose the 
correct 7, it is difficult to get a reasonable match between the experiment and theory. 



IV. CONTRIBUTION FROM LANDAU LEVEL 

In the Sondheimer oscillation we have pointed out that the predicted oscillation amplitude 
does not match with the experimental data because of the Shubnikov-de Haas oscillation 
which occurs above 4 T. In order to confirm irrelevance of the formation of Landau levels at 
low magnetic fields, we take into account the vector potential for the MR and Hall resistance. 
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FIG. 4: Hall resistance for a. jjl, b. [j2, and c jjS, where the same strength of disorder is utilized as 
in Fig. 3, respectively. We believe that the origin of this mismatch between experiment and theory 
lies in the contribution from the bulk transport which is not included in the single Dirac-fermion 
theory. 



This has been performed in the context of the quantum Hall effect in graphene, where the 
two valley contributions are simply added 



The longitudinal conductance is given by 
e'^NfV , 1 



X 



2u^ + 



UiUJ 



duo 



cosh^(^) (±^^)%(2c.r)^ 



A2 + r2 



r 



(a;2 _ A2 - r2)2 + 4a;2r2 
A2 + r2-u;2-2iu;r\i 
2v)\eH\/c J 



(6) 



with Nf = 1 {Nf = 2 for graphene), where A = hui is the Zeeman energy, F is the imaginary 
part of the electron self-energy due to disorder, and "^{z) is the digamma function. 
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FIG. 5: Each contribution in Hall conductance shows that the dominant term is the side-jump 
term in low fields and the Berry curvature term in high fields. These curves are generated, by 
using the parameter values from sample (13, a = 0.5 and 7 = 0.44. 




FIG. 6: a. Hall resistance deviates from the a = 0.7 line with 7 = 0.44, changing a. b. Hall 
resistance separates from the 7 = 0.44 line with a = 0.7, varying 7. The Hall resistance of sample 
[j4 is well fitted by using a = 0.7 and 7 = 0.44. 



The Hall conductance is 
where ub is the filling factor, given by 

Ub = 



e^Nf 



(7) 



du /CO + fi 
— tanh — — — 
27r V 2T 

T 

^V(a;-M„)2 + r2 

n=l ^ ' 



+ {u < — > -u) 



+ (cu < — > -cuj^ 



(8) 



where M„ = y + 2nf^|ei/|/c is the dispersion of surface Dirac electrons in the presence 
of the Landau level. 




FIG. 7: Quantized Hall conductance, expected to be relevant in high magnetic fields. 

By using these expressions, we plot the Hall conductance in Fig. 7. First of all, the 
plateau in the Hall conductance is clearly shown. This behavior is far from that given by 
the experimental data. Therefore, we conclude that our regime is far from being described 
by the quantum Hall effect. The introduction of the Zeeman term is sufficient to explain 
the transport data at low magnetic fields. 
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